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A  PAGE  TEST  WITH  NUISANCE  PARAMETER  ESTIMATION 

INTRODUCTION 

Signal  detectors  may  be  described  as  parametric,  where  the  observed  data  are  assumed  to 
follow  a  parameterized  probability  density  function  (PDF),  or  nonparametric.  Parametric 
detectors  provide  performance  improvement  over  nonparametric  detectors  at  the  expense  of 
increased  sensitivity  to  an  incorrect  parameterization.  Many  parameterizations  of  realistic 
detection  problems  result  in  composite  hypothesis  tests  where  there  are  unknown  nuisance 
and  signal  parameters.  Nuisance  parameters  are  unknown  parameters  having  the  same  value 
under  both  the  signal-present  and  signal-absent  hypotheses. 

When  testing  a  finite  set  of  data,  there  are  several  techniques  for  dealing  with  nuisance 
parameters:  uniformly  most  powerful  tests,  Bayesian  approaches  when  prior  distributions 
on  the  nuisance  parameters  are  known  or  can  be  assumed,  invariant  tests,  the  generalized 
likelihood  ratio  test,  or  ad  hoc  techniques  such  as  the  substitution  of  maximum  likelihood 
estimates  (MLEs)  of  the  nuisance  parameters  from  auxiliary  data  in  a  log-likelihood  ratio 
(LLR)  or  locally  optimal  nonlinearity.  Such  techniques  do  not  always  extend  to  sequential 
detection  problems.  As  illustrated  by  Govindarajulu  [1],  most  sequential  tests  for  composite 
hypotheses  essentially  estimate  unknown  parameters  using  all  data  previous  to  the  current 
sample.  The  disadvantages  of  these  techniques  are  that  they  do  not  allow  nonstationary 
nuisance  parameters  and  are  potentially  difficult  to  analyze.  Stahl  [2]  used  a  similar  pro¬ 
cedure  in  the  Page  test  [3]  for  a  change  detection  problem  where  the  nuisance  parameters 
were  estimated  by  an  exponential  averager  over  all  data  previous  to  the  current  sample. 
Disadvantages  of  this  technique  are  the  lack  of  theoretical  analysis  of  the  performance  and 
the  corruption  of  the  nuisance  parameter  estimate  (NPE)  by  data  containing  both  signal  and 
noise.  This  report  proposes  and  analyzes  a  method  for  the  estimation  of  nuisance  parameters 
in  the  Page  test  for  the  change  detection  problem. 
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BACKGROUND  -  THE  PAGE  TEST 


The  change  detection  problem,  in  general,  may  be  characterized  by  a  model  assuming  that 
the  observed  data 

Xi,X2,X3,...  (1) 


are  statistically  independent  and  distributed  according  to  the  rule 


Xi  ~  < 


/{x|A,^o) 

fix\X,0r) 


i  <p 
i>p 


(2) 


where  f{x\X,B)  is  a  parameterized  PDF  having  A  as  an  unknown  nuisance  parameter  and  9 
as  a  signal  strength  parameter.  Here,  is  the  value  of  $  when  no  signal  is  present  (typically 
$0  =  0),  9i  is  the  value  of  9  when  signal  is  present,  and  p  is  the  unknown  starting  time  of 
the  signal. 

The  change  detection  problem  amounts  to  determining,  as  rapidly  as  possible,  when  the 
sequence  of  data  changes  from  following  the  signal-absent  probability  distribution  law  to 
following  the  signal-present  law.  If  g  (x)  is  the  detector  nonlinearity  applied  to  the  observed 
data,  the  Page  test  declares  a  detection  when  the  Page  test  statistic. 


Wk  =  max{0,Wk-i+ g{xk)}  ,  k>l, 


(3) 


is  greater  than  a  threshold  h,  where  Wq  =  0.  A  sample  Page  test  statistic  sequence  is  found 
in  figure  1. 

As  with  most  sequential  tests,  the  primary  performance  measures  for  the  Page  test  are  the 
average  number  of  samples*  before  a  threshold  crossing  under  the  signal-present  and  signal- 
absent  hypotheses.  For  the  Page  test,  these  are  respectively  called  the  average  number  of 
samples  before  a  detection,  D,  and  the  average  number  of  samples  between  or  before  a  false 
alarm,  T. 

*Also  known  as  the  average  sample  number,  average  stopping  time,  or  average  record  length. 
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Figure  1.  Sample  Page  Test  Statistic  Sequence 


When  g  (x)  is  the  LLR  of  a  data  sample,  the  Page  test  yields  the  smallest  worst  case  D 
when  f  is  lower  bounded  by  Tq.  The  worst  case  D  is  observed  when  the  signal  starts  at  a 
regulation  or  reset  of  the  Page  test  statistic  (i.e.,  H4  =  0).  Had  the  signal  started  at  some 
other  time  index,  the  Page  test  statistic  may  have  crossed  the  threshold  sooner,  resulting  in 
a  smaller  value  for  D.  The  optimality  was  shown  asymptotically  as  Tq  oo  by  Lorden  [4] 
and  subsequently  for  finite  Tq  by  Moustakides  [5]. 

If  the  signal-present  hypothesis  is  characterized  by  an  unknown  signal  strength,  the  locally 
optimal  nonlinearity  may  be  substituted  for  the  LLR  in  g{x).  Dyson  [6]  showed  that  this 
results  in  the  same  asymptotic  optimality  (i.e.,  Lorden’s  result)  as  the  LLR  in  the  Page  test 
as  the  signal  strength  goes  to  zero. 

The  Page  test  has  also  been  applied  to  the  detection  of  a  temporary  change  in  the  PDF 
of  the  observed  data.  In  this  case,  the  false  alarm  performance  measure  (f)  remains  the 
same.  However,  the  detection  performance  is  characterized  by  the  probability  of  detecting 
the  signal.  Analytical  methods  for  approximating  the  probability  of  detection  that  are 
extremely  accurate  when  g  (xk)  is  Gaussian  are  discussed  by  Han,  Willett,  and  Abraham  [7]. 
A  simulation  based  method  providing  accurate  approximation  for  non-Gaussian  Page  test 
updates  may  be  found  in  [8]. 
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AVERAGE  SAMPLE  NUMBER  ANALYSIS 

Basseville  and  Nikiforov  [9]  discuss  the  primary  analytical  results  for  determining  the 
average  sample  numbers  (ASNs)  for  the  Page  test.  These  are  based  on  representing  the 
Page  test  as  consecutive  Wald  tests  terminating  at  the  signal-absent  or  null  hypothesis, 
followed  by  a  final  Wald  test  terminating  at  the  signal-present  hypothesis.  The  termination 
of  the  Wald  tests  at  the  null  hypothesis  is  reflected  in  the  Page  test  statistic  as  a  regulation 
or  reset  to  zero. 

A  Wald  test  may  be  described  by  the  stopping  time 

A  =  inf  {n  >  0  :  Sn  ^  (a,  b)}  (4) 


where  a  and  b  are  thresholds  chosen  to  satisfy  Type  I  and  II  error  probabilities  and 


Sn 


1=1 


(5) 


is  the  cumulative  summation  of  the  data  transformed  by  the  nonlinearity  g.  When 

g  (x)  is  the  LLR,  the  test  is  known  as  the  sequential  likelihood  ratio  test  (SLRT)  or  sequential 
probability  ratio  test  (SPRT)  and  is  optimal  in  the  sense  of  providing  the  smallest  ASNs  of 
any  test  with  equivalent  Type  I  and  II  error  probabilities  [10]. 

The  ASN  analysis  of  the  Page  test  employs  Wald’s  flrst  equation  ([11],  [9],  [12]) 


E[N]  = 


E[5^] 


(6) 


and  Wald’s  fundamental  identity 


E 


J.Sn 


M(t) 


-N 


=  1 


(7) 


where 

M(t)  =  E[e‘»(’’’^)]  (8) 

is  the  moment  generating  function  (MGF)  of  the  transformed  data.  Equation  (6)  is  used 
to  approximate  the  average  number  of  samples  for  the  Wald  test  terminating  at  either  the 
signal-present  or  signal-absent  hypotheses,  and  equation  (7)  is  used  to  approximate  the 
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probability  that  the  Wald  test  terminates  at  the  signal-absent  hypothesis  as  a  function  of 
the  signal  strength  6  (also  known  as  the  operating  characteristic  of  the  Wald  test).  How  the 
approximations  are  carried  out  forms  the  differences  between  the  techniques  discussed  by 
Basse ville  and  Nikiforov  [9]. 


Wald  Approximation 


The  first,  less  accurate  result  utilizes  Wald’s  approximations  to  the  ASNs  and  operating 
characteristic,  P  (d),  of  a  Wald  test.  They  are  based  on  the  assumption  that  when  the  test 
stops,  the  value  of  the  test  statistic  is  close  to  the  value  of  the  threshold;  that  is,  Sn  ^  o, 
or  Sn  ~  b.  As  seen  in  Basseville  and  Nikiforov  [9],  this  method  leads  to  the  approximations 
for  the  Page  test  ASNs 


and 


D 


1  +  hte^  - 
teo^eo  [g  (xi)] 


1  +  hU,  - 


(9) 


(10) 


te.Ee^igixi)] 

where  the  expectations  are  taken,  respectively,  under  the  signal-absent  (6*0)  and  signal-present 
(6>i)  hypotheses.  Here,  te,  for  9  =  60  or  9i,  is  the  nonzero  unity-root  of  the  MGF 


Eg  [e*''  =  M  {tg)  =  1  . 


(11) 


Siegmund  Approximation 

The  second  approximation  discussed  by  Basseville  and  Nikiforov  [9]  utilizes  corrections, 
derived  by  Siegmund  [11],  to  the  Wald  approximations  to  the  ASNs  and  operating  character¬ 
istic  of  the  Wald  test.  The  corrections  attempt  to  account  for  the  excess  above  or  below  the 
threshold  in  the  Wald  test  at  termination.  The  resulting  approximations  to  T  and  D  ([11], 
[9])  are  identical  in  form  to  those  of  equations  (9)  and  (10)  with  the  threshold  h  replaced  by 

hs  =  h  +  p+-p-  (12) 

where  p+  and  p_  represent  the  corrections  due  to  excess  beyond  a  threshold. 
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The  corrections  are  derived,  as  seen  in  Siegmund  [11],  by  first  normalizing  the  Wald  test 
statistic  to  have  zero-mean  and  unit-variance.  If  the  mean  and  variance  of  g  (x)  are  /i  and 
cr^  respectively,  then  the  normalized  update  is 

y.  =  ,  (13) 

a 

Let  the  cumulative  summation  be  described  by 

Vn  =  E^’ 

i=l 

and  let  the  test  (that  stops  when  the  upper  threshold  h  is  crossed)  have  the  stopping  time 

=  inf  (n  >  0  :  Vn'>  h}  .  (15) 


A  lemma  described  by  Siegmund  [13]  states  that  the  excess  beyond  the  threshold  in  the  Wald 
test  statistic  at  termination  may  be  approximated  by 


E 


In  [11],  Siegmund  shows  that 


where 


N+ 


E 

2E 

$  (A)  =  E 


E  [v)5,; 

(16) 

2E  [VnJ 

K  H 

(17) 

6  TT 

[l-$(A)]}dA, 

(18) 

'AVi] 

(19) 

is  the  characteristic  function  of  Yi, 

K  =  E[y*] 

=  (20) 

<7^ 

is  the  third  moment  of  Y  or,  equivalently,  the  skewness  of  g  (xi),  and  the  expectations  are 
taken  under  either  the  signal-present  or  signal-absent  hypotheses,  depending  on  whether  the 
correction  is  for  D  or  T. 
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The  excess  beyond  the  lower  threshold  will  simply  be  the  negative  of  equation  (17)  if  the 
Wald  test  update  has  a  symmetric  distribution  as  is  the  case  with  a  Gaussian  g  (xi)  (this 
is  all  that  Siegmund  considers  in  [11]).  Let  the  test  that  stops  when  the  lower  threshold  is 
crossed  have  the  stopping  time 

iV_  =  inf{n>0:  W  <  0}  .  (21) 


The  excess  in  the  Wald  test  statistic  at  termination  may  then  be  approximated  by 


where 


E 


N- 


e[W| 


2E  [l 

E 

2' 

2E  [-Vn_ 

E 

K4 

2E 

Wm^[ 

W„  =  -V„  =  J2i-Y’) 

i=l 


and 


(22) 


(23) 


M+  =  inf  {n  >  0  :  Wn  >  0}  =  iV_  , 


(24) 


which  is  the  same  form  as  the  above  (equations  (14)-(16))  except  that  Yi  is  replaced  by 
—Yi-  Noting  that  the  third  central  moment  of  —Yi  is  simply  —n  and  that  the  characteristic 
function  of  -Yi  is  $  (-A),  equation  (22)  may  be  simplified,  using  equation  (17),  to 


2E 

Wm^' 

^  ‘os  (I  [1  -  » • 


(25) 


The  natural  logarithm  of  a  complex  number  {z  =  reP^)  may  be  described  in  Cartesian 
form  through  the  following  identity 


log  z  =  log  re^® 

=  logr  +  j^. 


(26) 
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Thus,  the  real  part  of  the  natural  logarithm  of  a  complex  number  is  simply  the  logarithm  of 
the  magnitude, 


Relogz  =  logr.  (27) 

Use  of  this  identity  in  equation  (18)  results  in  the  simplification 

=  A-Uog{^|l-*(A)l}(iA.  (28) 

Now,  using  the  identity 


$(A)  =  E 


J^Vi 


=  E  [cos  (AT))]  +  jE  [sin  (ATj)]  , 


it  can  be  shown  that 

|1-$(A)|  =  |1  -  E[cos(Ayi)] -iE[sin(Ari)]| 

=  {(1  -  E  [cos(Ay;)])^  +  E  (sin  (Ay.)]"}' 

=  {(1  -  E [cos (-Ay)])"  +  E [sin (-Ay)]"}' 
=  |l-<!>(-A)| 


(29) 


(30) 


is  an  even  function  of  A.  Thus,  equation  (25)  simplifies  to 

_E[l^  ^  Py 

2E  [WmJ  6  TT 

Incorporating  a  non-unit-variance  into  the  Wald  test  update  simply  results  in  scaling 
equations  (17)  and  (31)  by  the  standard  deviation  a  to  yield  the  Siegmund  corrections 


(32) 


and 
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EXAMPLES 

Gaussian  shift  in  mean,  exponential,  and  noncentral  chi-squared  signal  types  often  arise 
in  sonar  signal  processing  applications.  For  example,  the  Gaussian  shift  in  mean  signal  may 
represent  a  range-induced  bias  in  the  innovations  sequence  of  a  tracking  algorithm.  The 
exponential  signal  type  is  encountered  as  the  magnitude-squared-basebanded-matched-filter 
output  or  as  the  result  of  frequency  domain  processing  on  a  sinusoidal  signal  with  Gaussian 
random  amplitude  and  uniform  random  phase  corrupted  by  Gaussian  noise  (e.g.,  multiple 
point  reflectors  with  none  dominant  for  active  sonar).  The  noncentral  chi-squared  signal  is 
the  result  of  frequency  domain  processing  on  a  sinusoidal  signal  with  deterministic  amplitude 
and  uniform  random  phase  corrupted  by  Gaussian  noise  (e.g.,  a  single  dominant  reflector 
for  active  sonar).  In  the  following  sections,  the  Siegmund  correction  terms  for  each  of  these 
signal  types  are  presented  and  used  to  compare  the  Siegmund  and  Wald  approximations  to 
D  and  f  to  the  ASN  obtained  through  simulation. 

Siegmund  Correction  Terms 

This  section  contains  the  Siegmund  correction  terms  and  the  information  required  to 
obtain  them  for  Gaussian  shift  in  mean,  exponential,  and  noncentral  chi-squared  signals. 
Table  1  contains  the  Siegmund  correction  terms  and  the  mean,  variance,  skewness,  and 
characteristic  function  of  the  normalized  statistic  of  equation  (13)  for  each  of  the  assumed 
distributions  for  the  Page  test  update  g{x).  Here,  represents  a  Gaussian  random 

variable  with  mean  /i  and  variance  Exp  (a)  represents  an  exponential  random  variable 
with  mean  a,  and  {5)  represents  a  noncentral  chi-squared  random  variable  with  n  degrees 
of  freedom  and  noncentrality  parameter  6. 

The  results  for  the  Gaussian  shift  in  mean  signal  may  be  found  in  Siegmund  [11]  or 
Basseville  and  Nikiforov  [9].  The  skewness  values  and  the  characteristic  functions  for  the 
unnormalized  random  variables  may  be  found  in  Manoukian  [12]  or  Johnson  and  Kotz  [14]. 
The  Siegmund  corrections  for  the  noncentral  chi-squared  signal,  p*^  and  p*_,  are  functions  of 
n  and  5  and  must  be  computed  individually  using  equations  (32)  and  (33).  Table  2  contains 
numerical  values  of  the  correction  terms  for  n  =  2  and  several  values  of  6.  Note  that  when 
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n  —  2  and  <5  =  0,  a  central  chi-squared  random  variable  results  and  the  correction  terms 
are  identical  to  those  for  the  exponential  signal  with  cr  =  2,  which  correctly  reflects  the 
relationship  between  these  two  types  of  random  variables. 


Table  1.  Siegmund  Corrections  to  Sample  Signals 


EM 

Var  [g] 

K 

$(A) 

P+ 

P- 

0 

a2 

e  2 

0.5825a 

-0.5825a 

Exp {a) 

a 

2 

e-tA(i  -  j\)~^ 

a 

{5) 

n-\-  6 

0-2  =  2  (n  -I-  2^) 

\/8{n+35) 

P\ 

o* 

(n+25)i 

Table  2.  Siegmund  Corrections  for  Noncentral  Chi-Squared  Signal  for  n  =  2 

and  Several  Values  of  6 


s 

Pi 

P- 

0 

2 

2 

3 

0.1 

2.0948 

-0.7008 

1 

2.6528 

-0.9861 

10 

4.9129 

-2.9735 

Simulation  Comparison 

Ideally,  the  nonlinearity  applied  to  the  Page  test  is  the  LLR  or  the  locally  optimal  non¬ 
linearity  of  the  observed  data.  For  the  Gaussian  shift  in  mean  and  exponential  signal  types, 
the  LLR  and  locally  optimal  nonlinearity  have  the  same  form:  scaled  shifted  versions  of  the 
data.  The  LLR  for  the  noncentral  chi-squared  signal  has  a  complicated  form  [14];  however, 
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the  locally  optimal  nonlinearity  is  simply  the  shifted  data.  The  amount  of  the  shift  may  be 
chosen  to  maximize  the  asymptotic  performance  of  the  Page  test  for  a  specified  signal-to- 
noise  ratio  (SNR)  as  described  in  [15].  Table  3  contains  the  statistical  description,  LLR  or 
locally  optimal  nonlinearity,  bias  term,  mean,  variance,  and  SNR  definition  for  each  of  the 
signal  types. 


Table  3.  Information  About  Sample  Signals 


Type 

9{x) 

Bias  (r) 

E.bl 

Var^  [p] 

SNR 

H 

1 

U 

xL 

2(t2 

id 

20-2 

d 

d 

Exp  (1  -t-  s) 

log  (1  +  s) 

S  —  T 

s 

w 

X  —  T 

2  (l  -b  1)  log  (l  -b  f ) 

6  +  2-t 

4((5-bl) 

6 

In  order  to  determine  the  Wald  and  Siegmund  approximations  to  D  and  T,  the  unity-root 
of  the  MGF  (11)  is  required  under  the  signal-present  and  signal-absent  hypotheses.  As  the 
nonlinearity  for  the  Gaussian  shift  in  mean  and  exponential  signals  is  their  respective  LLR, 
the  MGF  roots  are  too  =  1  and  ^  -1.  The  MGF  roots  for  the  locally  optimal  nonlinearity 
for  the  noncentral  chi-squared  signal  are  the  solutions  to 

- Ttei  -  log  (1  -  2fei)  0  (34) 

and 

-b  log  (1  -  2%)  =  0  .  (35) 

The  Wald  and  Siegmund  approximations  for  the  average  delay  before  detection  and  that 
estimated  from  simulations  are  found  in  figures  2-4  for  each  signal  type  as  a  function  of 
SNR.  The  threshold  was  /i  =  5  for  the  Gaussian  shift  in  mean  and  exponential  signals 
and  h  =  10  for  the  noncentral  chi-squared  signal.  The  estimated  value  of  D  is  the  sample 
mean  of  the  number  of  samples  required  for  detection,  taken  over  1000  trials.  Observe 
that  for  each  signal  type,  the  Siegmund  approximation  is  closer  to  the  simulated  value  than 
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the  Wald  approximation,  particularly  for  large  SNR.  This  result  is  expected,  as  the  Wald 
approximation  ignores  the  excess  over  the  threshold,  which  may  be  substantial  for  large 
strength  signals. 

The  Wald  and  Siegmund  approximations  for  the  average  time  between  false  alarms  and 
that  estimated  from  simulations  are  found  in  figures  5-7  for  each  signal  type  as  a  function  of 
threshold.  The  nonlinearity  for  each  signal  was  formed  assuming  that  the  signal  strength  was 
zero  decibels  (^  =  s  =  5  =  l).  As  with  the  results  for  D,  utilizing  the  Siegmund  correction 
terms  provides  a  better  approximation  than  the  straight  Wald  method.  The  quality  of  the 
Wald  and  Siegmund  approximations  is  not  as  good  for  T  as  for  D.  This  is  a  direct  result  of 
the  increased  number  of  resets  or  regulations  of  the  Page  test  when  no  signal  is  present  (i.e., 
more  Wald  tests  terminating  at  the  signal-absent  hypothesis),  which  introduces  a  higher 
sensitivity  to  errors  in  the  approximations  to  the  ASNs  and  operating  characteristic  of  the 
Wald  test. 
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Figure  2.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Delay  Before  Detection  for  Gaussian  Shift  in  Mean  Signal  With  h  =  5 


Figure  3.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Delay  Before  Detection  for  Exponential  Signal  With  /i  =  5 
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Figure  4.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Delay  Before  Detection  for  Noncentral  Chi-Squared  Signal  With  h=  10 


Figure  5.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Time  Between  False  Alarms  for  Gaussian  Shift  in  Mean  Signal 
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Figure  6.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Time  Between  False  Alarms  for  Exponential  Signal 


Figure  7.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Time  Between  False  Alarms  for  Noncentral  Chi-Squared  Signal 
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PROPOSED  ALGORITHM 

As  previously  mentioned,  the  Page  test  [3]  may  be  described  as  consecutive  Wald  tests  ter¬ 
minating  at  the  signal-absent  hypothesis  followed  by  a  final  Wald  test  terminating  at  the 
signal-present  hypothesis.  Each  termination  of  a  Wald  test  at  the  signal-absent  hypothesis 
(represented  by  a  regulation  or  reset  of  the  Page  test  to  zero)  indicates  that  all  previous 
data  are  considered  to  be  signal  free.  Thus,  these  data  may  be  used  to  estimate  any  un¬ 
known  parameters  that  also  exist  under  the  signal-absent  hypothesis;  for  example,  nuisance 
parameters.  An  additional  benefit  of  using  the  data  previous  to  the  most  recent  reset  of  the 
Page  test  is  the  statistical  independence  of  these  data  and  the  Page  test  statistic  after  the 
reset.  Thus,  analysis  may  be  performed  by  conditioning  on  the  observed  NPE,  followed  by 
an  expectation  removing  the  conditioning. 

The  following  update  algorithm  describes  the  proposed  Page  test  with  nuisance  parameter 
estimation.  Here,  auxiliary  data  are  data  that  have  been  determined  to  be  signal-free  by 
the  Page  test,  A  is  the  NPE  formed  from  the  auxiliary  data,  and  (x)  is  the  detector 
nonlinearity. 

1.  If  Wfc  =  0 

-  Store  new  auxiliary  data 

-  Update  A  using  new  auxiliary  data 

2.  Set  Wfc+i  =max{0,H4-|-^^(xfc+i)} 

3.  If  Wk+i  >  h 

-  Declare  signal  present  and  stop 
Else 

-  Set  k  =  k  +  1  and  return  to  step  1 

In  the  following  sections,  details  associated  with  the  nuisance  parameter  estimation  are 
discussed  and  the  Siegmund-  and  Wald-based  approximations  to  the  ASN  of  the  Page  test 
with  nuisance  parameter  estimation  are  derived  and  shown  to  simplify  to  those  for  the  stan¬ 
dard  Page  test  when  the  nuisance  parameter  is  known  exactly.  The  Page  test  with  nuisance 
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parameter  estimation  is  then  applied  to  detect  a  Gaussian  shift  in  mean  signal  where  D 
and  T  are  derived  and  compared  to  simulation  results.  The  Siegmund-based  approximations 
are  used  to  illustrate  the  Page  test  performance  as  a  function  of  the  threshold,  SNR,  and 
amount  of  auxiliary  data.  Finally,  the  Wald-based  approximations  are  used  to  illustrate  the 
asymptotic  (in  the  sense  of  large  threshold)  relationship  between  the  threshold  and  f  and 
D. 

NUISANCE  PARAMETER  ESTIMATION 

The  NPE  is  formed  from  auxiliary  data  and,  depending  on  the  situation,  may  only  be 
a  consistent  estimator  under  the  signal-absent  hypothesis.  Thus,  it  is  important  that  the 
auxiliary  data  truly  be  signal-free.  A  necessary  but  not  sufficient  condition  for  ensuring  this 
is  to  have  the  auxiliary  data  consist  of  data  previous  to  the  most  recent  reset  of  the  Page 
test.  However,  due  to  the  inherent  delay  observed  between  the  time  a  signal  starts  and  a 
threshold  crossing  indicating  a  detection,  there  is  a  nonzero  probability  that  this  auxiliary 
data  may  contain  both  signal  and  noise.  To  counter  this  possibility,  it  is  suggested  that  a 
fixed  number  of  either  samples  or  resets  be  used  as  a  buffer  zone.  Figure  8  depicts  these  two 
arrangements. 

For  extremely  weak  signals,  the  delay  before  detection  can  be  substantial.  This  may  cause 
the  NPE  to  be  corrupted  by  data  containing  noise  and  signal  even  with  a  carefully  designed 
buffer  zone.  This  problem  may  be  avoided  altogether  if  there  exists  an  estimator  for  the 
nuisance  parameter  that  has  the  same  form  under  both  the  signal-absent  and  signal-present 
hypotheses. 

Most  practical  scenarios  will  involve  nuisance  parameters  that  are  slowly  changing  in 
time.  In  this  situation,  the  auxiliary  data  should  consist  of  a  block  of  data — as  opposed  to 
using  all  previous  data  to  estimate  the  nuisance  parameter.  The  size  of  the  block  and  its 
positioning  should  be  chosen  to  ensure  that  the  nuisance  parameter  is  stationary  from  the 
beginning  of  the  block  of  auxiliary  data  to  either  a  threshold  crossing  indicating  a  detection 
or  the  next  reset  of  the  Page  test  statistic. 
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Page  Test 

Auxiliary  Buffer  /  Statistic 

Data  Zone  "  •  -  ■' 


m  Samples  Before  Most  Recent 
Most  Recent  Reset  Reset 


Page  Test 

Auxiliary  Buffer  ,  .•  Statistic 

Data  Zone 


mth  Most  Recent 

Recent  Reset  Reset 


Figure  8.  Placement  of  Auxiliary  Data  in  Relation  to  the  Most  Recent  Reset  of 

the  Page  Test  Statistic 
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AVERAGE  SAMPLE  NUMBER  ANALYSIS 

As  seen  in  Basseville  and  Nikiforov  [9],  the  first  step  in  determining  the  ASN  for  the 
Page  test  is  to  determine  the  relationship  to  the  ASNs  and  operating  characteristic  for  the 
individual  Wald  tests  that  form  the  Page  test.  The  following  definitions  are  required; 


A  (9)  -  Average  number  of  samples  for  the  Page  test  with  nuisance  parameter  es¬ 
timation  when  9  is  the  signal  strength 

Bh  (9)  -  Average  number  of  samples  for  the  Wald  test  terminating  at  threshold  h 
when  9  is  the  signal  strength 

Bq  (9)  -  Average  number  of  samples  for  the  Wald  test  terminating  at  threshold  0 
when  9  is  the  signal  strength 

P  (9)  -  Probability  that  the  Wald  test  terminates  at  threshold  0  when  9  is  the 
signal  strength 

Suppose  the  number  of  Wald  tests  that  comprise  the  Page  test  is  m;  that  is,  m  -  1  Wald 
tests  terminating  at  threshold  0  followed  by  one  terminating  at  threshold  h.  Then,  the  ASN 
for  the  Page  test  conditioned  on  m  is 

A{9\m)  =  {m-1)  Bo  {9)  +  Bn  {9)  .  (36) 

The  conditioning  on  m  may  be  removed  by  taking  the  expected  value 

A{9)  =  E[A{9\m)] 

=  E[m-l]Bo{9)  +  Bn{9) 

=  E[m]Bo{9)  +  Bn{9)~Bo{9)  .  (37) 

Each  consecutive  Wald  test  in  the  Page  test  is  independent  of  the  others.  Thus,  as  found  in 
Manoukian  [12],  m  follows  a  geometric  probability  mass  function  with  parameter  q  =  1-P  {9) 
and  mean 


E  [m 


1 

q 

1 

1-F(«)  ■ 


(38) 
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This  results  in  the  form 


(39) 


for  the  ASN  of  the  Page  test. 

When  the  Page  test  includes  nuisance  parameter  estimation,  with  the  update  scheme 
described  in  the  previous  section,  the  ASNs  and  operating  characteristic  of  the  Wald  test 
must  be  determined  conditioned  on  the  NPE.  Here,  it  is  assumed  that  each  NPE  is  inde¬ 
pendent  of  the  data  in  the  nearby  Wald  tests  and  independent  of  all  the  other  NPEs.  This 
assumption  will  be  valid  at  least  for  large  buffer  sizes  and  short  block  sizes.  Now,  the  ASNs 
and  operating  characteristic  of  the  Wald  test  may  be  described  as 


bu«)  =  e4b»(9|a)]  . 


B„(«)  =  Ej[Bo(«]A)' 


(40) 

(41) 


and 

P(9)  =  Bs[b(9|a)]  .  (42) 

where  the  arguments  of  the  expectations  are  the  respective  ASN  or  operating  characteristic 
conditioned  on  the  NPE  A. 

The  test  statistic  of  the  Wald  test  used  to  implement  the  Page  test  at  the  sample  is 
the  cumulative  summation 

Sn  =  f^gx{^i)  ■  (43) 

Let  N  be  the  stopping  time  of  this  particular  Wald  test  (i.e.,  the  test  has  an  upper  threshold 
of  h  and  a  lower  threshold  of  zero). 


A  =  inf  {n  >  0  :  ^  (0,  h)}  . 


(44) 


As  found  in  Siegmund  [11]  or  [13],  the  average  excess  of  the  test  statistic  beyond  the  thresh¬ 
olds  at  termination  may  be  approximated  by 


h 


Sn  ^  h,X 


~  P+ 


(45) 
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and 


E 


S'/v  ‘S'iv  ^  0,  A 


P-  , 


(46) 


where  p+  and  p_  are  determined  in  equations  (32)  and  (33)  and  are  functions  of  the  NPE  A. 

Then,  using  Wald’s  first  equation  (6),  the  ASNs  of  the  Wald  test,  conditioned  on  the 
NPE  A,  may  be  written  as 

Bh{9\x)  -- 


E 

Sn 

Si\f  ^  h,X 

E 

axi^) 

X 

h  +  E  -  h ' 

Sn 

>  h,X 

E 

ax  (^) 

X 

h  +  P+ 


E 


[ax  (^) 


A 


(47) 


and 


b„(«|a) 


E 

Sn 

Sn  <  0,A 

E 

P 

axi^) 

X 

E 

ax  (^)  ^ 

(48) 


Following  Siegmund  [11]  or  Basseville  and  Nikiforov  [9],  the  operating  characteristic  of  the 
Wald  test  may  similarly  be  related  to  p+  and  p_,  except  that  here  it  is  conditioned  on  the 
NPE  A.  Let  the  conditional  MGF  of  the  detector  nonlinearity  g-^  (x)  be 

[e*  3a  W 


A 


(49) 


Under  very  mild  restrictions  on  (x),  it  can  be  shown  that  there  exists  a  nonzero  value  te 
such  that  the  MGF  is  one,  M-^  {te)  =  1,  where  0  is  the  signal  strength  parameter.  Substitution 
of  this  value  into  Wald’s  fundamental  identity, 


E 


N 


A 


=  1 


(50) 


results  in 


1  =  E 


Sn 


'i-p[e 

'  E 

Q^e  Sn 

SN>h,x]+P  {9 

a)e 

Sn 

Sn  ^  0,  a 

(61) 
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Solving  for  1  -  P  A)  results 


in 


p{e 


A  - 


1  -E 

e*"  Sn<0,X 

E  e‘9  Sn  >h,X 

1-E 

-e[< 

0*9  Sn 

^tg  Sn  ^ 

Sn  < 

Sn  <  0,A 

0,A 

6^  E  6^® 

Sn 

>  h,X 

,]  -E 

Qtg  Sn  Sn  <0,X 

TT  •  (52) 


Siegmund  [11]  and  Basseville  and  Nikiforov  [9],  in  analyzing  the  standard  Page  test,  now 
proceed  by  liberally  applying  a  first  order  approximation  to  the  exponential  function.  The 
same  results  may  be  obtained  by  recognizing  that  the  expectations  in  the  numerator  and 
denominator  of  equation  (52)  are  describable  as  the  MGF  of  a  random  variable  evaluated 
at  a  small  argument  {-Sn  is  expected  to  be  small  when  the  test  terminates  by  crossing  the 
threshold  at  zero  and  -  his  expected  to  be  small  when  the  test  terminates  by  crossing 
the  threshold  at  h).  Thus,  the  following  first  order  approximation  may  be  applied, 

E 


=  E 

[l+ar+<“7N...l 

2! 

^ ,,,,  a^E  [r^j 
=  1  +  oE  [y]  + - — -  + 


2! 


1  +  aE  [Y] 

eaElv] 


(53) 


Disregarding  terms  on  the  order  of  {tep-)  ,  equation  (52)  becomes 

i-p(e\x)  ? 
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[^t^E  Sn 

Sn  ^  0)  A 

1 

Qhtg  gxp  j 

Sn  h 

Al 

}  -  exp  j 

[ieY 

Sn 

5v<0,a]} 

1 
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-  1 
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_ (54) 

1  _  Qtg{h+p+-p-) 

Substituting  equations  (47),  (48),  and  (54)  into  equations  (40),  (41),  and  (42),  respectively, 
followed  by  substitution  into  equation  (39),  results  in  the  Siegmund-based  approximation  to 
the  ASN  of  the  Page  test  with  nuisance  parameter  estimation. 
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(55) 
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The  corresponding  Wald-based  approximation  ignores  the  excess  over  the  boundary;  that  is, 
->  0  and  p-  — >  0.  Applying  this  to  equation  (55)  and  utilizing  L’Hospital’s  rule  yields 
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(56) 


Thus,  equations  (55)  and  (56)  describe  the  Siegmund-  and  Wald-based  approximations  to 
the  ASN  of  the  Page  test  with  nuisance  parameter  estimation.  The  MGF  root  te,  Siegmund 
corrections  p+  and  p_,  and  nonlinearity  mean  E[gj^(x)]  must  be  determined  on  a  case-by¬ 
case  basis  as  functions  of  the  NPE  A  and  substituted  into  these  equations.  Under  special 
circumstances  (e.g.,  te^  <  0  and  teo  >  0  for  all  A)  these  equations  may  yield  simpler  forms 
accurate  for  approximation  to  T  and  D. 


Relationship  to  Standard  Page  Test 


If  the  amount  of  data  used  to  estimate  the  nuisance  parameter  A  tends  to  infinity  and 
the  NPE  A  is  consistent,  A  ^  A  and  the  Page  test  with  nuisance  parameter  estimation 
becomes  the  standard  Page  test.  Thus,  the  expectations  in  equations  (55)  and  (56)  over  A 
are  removed,  which,  as  expected,  results  in  the  standard  Siegmund-based  approximation 


^5  {9) 


l  +  to{h  +  p+  -  pJ)  — 
tpE  [gx  (x)] 


(57) 


and  the  standard  Wald-based  approximation 


Aw  (9) 


1  +  teh  —  e*®^ 
tpE  [gx  (x)] 


(58) 
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GAUSSIAN  SHIFT  IN  MEAN  SIGNAL  WITH  UNKNOWN  VARIANCE 


As  previously  mentioned,  a  signal  type  commonly  encountered  in  sonar  signal  processing 
may  be  modeled  by  a  change  in  the  mean  of  a  Gaussian  random  variable.  Additionally, 
suppose  that  the  variance  of  the  Gaussian  random  variable  is  unknown  but  constant.  Under 
these  conditions,  the  observed  data  are  distributed  according  to 


Xi  ~  < 


(59) 


A/’(0,(t^)  i<p 
A/'(/i,cr^)  i>p 

where  p  is  the  unknown  starting  time  of  the  signal,  the  mean  /i  is  known,  and  the  variance 
(7^  is  unknown  and  considered  to  be  a  nuisance  parameter.  If  the  auxiliary  data,  say  {Yj}j_^, 
are  assumed  to  be  signal-free,  the  MLE  of  cr^  is 

=  (60) 
In  this  case,  the  unbiased  estimator  for  the  variance  under  the  signal-present  hypothesis. 


1  •!+!  ,  ,2 


(61) 


may  be  equivalently  used,  where 


Y  = 


7=1 


1  J+1 


(62) 


and  it  is  assumed  that  the  estimate  is  formed  from  J  +  1  samples  of  auxiliary  data  rather 
than  J  to  retain  statistical  equivalence  between  the  estimators. 

Appropriately  scaled,  (for  either  of  the  above  estimators)  follows  a  central  chi-squared 
distribution  with  J  degrees  of  freedom, 


w  =  ^~xj. 

/T^ 


(63) 


Also  note  that  is  unbiased  for  that  is,  E  [d^]  =  cr^. 

The  LLR  or  locally  optimal  nonlinearity,  with  the  NPE  d^  substituted  for  has  the 

form 


9(^)  =  |^(^-|)  . 


(64) 
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2 

which,  in  this  case,  is  the  actual  LLR  scaled  by  The  mean,  variance,  and  MGF  unity-root 
of  g  (x)  conditioned  on  the  NPE  under  the  signal- present  (0  =  E  [x]  =  g)  and  signal-absent 
(^  =  E  [x]  =  0)  hypotheses  are  found  in  table  4.  Note  that  the  MGF  unity-roots  are  simply 
±1  (the  MGF  unity-roots  for  the  LLR)  divided  by  the  scale  term, 


Table  4.  Mean,  Variance,  and  MGF  Unity-Root  for  Gaussian  Shift  in  Mean 

Signal  Conditioned  on  the  NPE 


Signal-present 

Signal-absent 

{6  =  p) 

(9  =  0) 

Ee  [g  (x)] 

jI. 

2^2 

Ziit 

2^2 

Yave  [g  (x)] 

MGF  Root  te 

0-2 

(7^ 

Conditioned  on  d^,  5'  (x)  is  still  Gaussian  under  both  the  signal-present  and  signal-absent 
hypotheses.  Thus,  as  found  in  table  1,  the  Siegmund  correction  terms  are 

-  0.5825^  (65) 

and 

p-  =  -p+  =  -0.5825^  .  (66) 

Average  Delay  Before  Detection 

In  order  to  determine  the  Siegmund-based  average  delay  before  detection,  t)s  — 
As  {9i  =  p),  each  of  the  terms  in  equation  (55)  must  be  determined  under  the  signal-present 
hypothesis.  For  notational  convenience,  let  7  =  0.5825.  These  terms,  which  are  expectations 
over  the  NPE  (in  this  case,  d^),  are  as  follows: 


'h  +  p+  -p- 

=  E^2 

2d^  f,  ^  aG\ 
— —  (  h  -t-  27—  ) 

_E[p(x)  |d2]_ 

p2  V  G^J_ 
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(68) 


(69) 


where 


eK  =  E 


e-2fc7f 


(70) 


fc=i  (1  +  2A;5)  " 

is  used  to  approximate  the  infinite  summation  in  the  next  to  last  line  of  equation  (69).  To 
arrive  at  the  result  of  equation  (69),  the  MGF  of  the  chi-squared  random  variable  W, 
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w 


AW 


=  (1  -  2ty^ 


(71) 


is  required. 
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Substitution  of  equations  (67)-(69)  into  equation  (55)  results  in  the  Siegmund-based 
approximation  to  D  for  the  Gaussian  shift  in  mean  signal  with  unknown  variance, 


Ds  =  [ha‘^  +  2'y jia^  - 


IJ?  (1  +  Ck) 


Jii 

2<t2 


(72) 


In  order  to  determine  the  Wald-based  approximation  to  the  average  delay  before  detec¬ 
tion,  Dw  =  Aw  {0i  =  Z^),  the  following  terms  in  equation  (56)  must  be  determined  under 
the  signal-present  hypothesis: 


and 
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(74) 


where 
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K 

<5k  =  E- 

fe=i  (l  -I-  2k^) 


:  +  l 


(75) 


is  used  to  approximate  the  infinite  summation  in  the  next  to  last  line  of  equation  (74). 
To  arrive  at  the  result  of  equation  (74),  the  first  derivative  of  the  MGF  of  the  chi-squared 
random  variable  W  is  required, 
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-  J(l-2i)'2"V  (76) 


Substitution  of  equations  (73)  and  (74)  into  equation  (56)  results  in  the  Wald-based 
approximation  to  D  for  the  Gaussian  shift  in  mean  signal  with  unknown  variance, 


Dw  — 


h- 


1+Sk 


(77) 


The  results  of  equations  (72)  and  (77)  are  compared  to  the  simulated  average  delay  before 
detection  as  a  function  of  SNR  in  figure  9  for  J  =  20  and  h  =  5.  The  simulation  used  1000 
trials  that  were  terminated  if  a  threshold  crossing  did  not  occur  within  1000  samples  (all 
trials  terminated  prior  to  this  limit).  Throughout  this  and  the  following  sections,  the  6^  and 
6k  summations  in  the  Siegmund  and  Wald  approximations  are  terminated  at  K  =  10  and 
the  variance  estimator  of  equation  (61)  is  used.  Similar  to  the  analysis  of  the  Page  test  when 
the  nuisance  parameters  are  known,  the  Siegmund-based  approximation  of  equation  (72) 
provided  substantially  more  accurate  results  than  the  Wald-based  approximation. 


Average  Time  Between  False  Alarms 

In  order  to  determine  the  Siegmund-based  average  time  between  false  alarms,  Ts  — 
As  {Bq  =  0),  each  term  in  equation  (55)  must  be  determined  under  the  signal-absent  hypoth¬ 
esis: 
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and 
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Substitution  of  equations  (78)-(80)  into  equation  (55)  results  in  the  Siegmund-based  approx¬ 
imation  to  f  for  the  Gaussian  shift  in  mean  signal  with  unknown  variance, 


T5  =  -^  (ha‘^  +  27/icr)  -h 


_j£i 

2(72 


(81) 


In  order  to  determine  the  Wald-based  approximation  to  the  average  time  between  false 
alarms,  fw  =  Aw  {Oq  =  0),  each  of  the  terms  in  equation  (56)  must  be  determined  under 
the  signal-absent  hypothesis: 
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1  OO 

k=i 


-Sk 


(83) 


Substitution  of  equations  (82)  and  (83)  into  equation  (56)  results  in  the  Wb  Id-based  approx¬ 
imation  to  T  for  the  Gaussian  shift  in  mean  signal  with  unknown  variance, 

fw  =  (84) 

_  rl- 

2^2 

The  results  of  equations  (81)  and  (84)  are  compared  to  the  simulated  average  time  be¬ 
tween  false  alarms  as  a  function  of  threshold  in  figure  10  for  J  =  20  and  a  zero  dB  SNR. 
The  simulation  used  5000  trials  that  were  terminated  if  a  threshold  crossing  did  not  occur 
within  2500  samples  (two  trials  failed  to  terminate  prior  to  this  limit).  Again,  the  Siegmund- 
based  approximation  of  equation  (81)  provided  substantially  more  accurate  results  than  the 
Wald-based  approximation. 


Page  Test  Performance 

The  previous  two  sections  have  demonstrated  the  quality  of  the  Siegmund-based  approx¬ 
imations  to  D  and  T  for  the  Page  test  with  nuisance  parameter  estimation.  In  this  section, 
these  approximations  are  used  to  investigate  the  performance  as  the  amount  of  auxiliary 
data,  the  SNR,  and  the  threshold  are  varied.  Figure  11  contains  a  Page  test  receiver  operat¬ 
ing  characteristic  (a  plot  of  D  against  f)  where  it  is  seen  that  the  performance  tends  to  that 
of  the  standard  Page  test  when  the  variance  of  the  data  is  known  ( J  =  oo)  as  the  amount 
of  auxiliary  data  increases.  Figure  12  demonstrates  how  the  threshold  required  to  achieve 
a  desired  f  increases  as  the  amount  of  auxiliary  data  diminishes,  which  is  explained  by  the 
increased  variability  of  the  Page  test  update  (xj)  due  to  the  increased  variability  of  the 
NPE  A.  The  loss  incurred  by  having  to  estimate  the  unknown  variance  is  seen  in  figure  13 
as  an  increase  in  D  when  T  is  held  constant. 
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Figure  9.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Delay  Before  Detection  for  Gaussian  Shift  in  Mean  Signal  With  Threshold 
=  5  for  Page  Test  With  Nuisance  Parameter  Estimation 


Figure  10.  Simulation  and  Wald  and  Siegmund  Approximations  to  the  Average 
Time  Between  False  Alarms  for  Gaussian  Shift  in  Mean  Signal  for  Page  Test 

With  Nuisance  Parameter  Estimation 
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Figure  11.  Average  Delay  Before  Detection  Vs.  Average  Time  Between  False 
Alarms  for  0  dB  SNR  and  Varying  Amounts  of  Auxiliary  Data 


Figure  12.  Average  Time  Between  False  Alarms  Vs.  Threshold  for  0  dB  SNR 

and  Varying  Amounts  of  Auxiliary  Data 
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Asymptotic  Performance 


Most  applications  of  the  Page  test  will  require  a  large  average  time  between  false  alarms. 
This  results  in  a  large  threshold,  h,  which  can  simplify  the  computation  of  the  threshold 
required  to  achieve  a  specific  f  and  also  provide  insight  into  the  relationship  between  D  and 
f.  Assuming  that  h  is  large,  the  Wald  approximations  to  f  and  D  for  the  standard  Page 
test,  equations  (9)  and  (10),  result  in  the  asymptotic  approximations 


T 


-tOo^Oo  [9  (2:)] 


and 


D 


h 


(85) 


(86) 


[9  (2:)]  ’ 

where  the  conditions  >  0  and  <  0  have  been  exploited.  The  exponential  relationship 
between  h  and  f  and  the  linear  relationship  between  h  and  D  are  evident  from  equations  (85) 
and  (86). 

Broder  [16]  carried  this  analysis  a  step  further  to  demonstrate  the  asymptotic  (large 
h)  exponential  relationship  between  f  and  D  by  introducing  an  asymptotic  performance 
measure  and  relating  it  to  the  MGF  root  under  the  signal-absent  hypothesis  and  the  mean 
of  the  Page  test  statistic  update  under  the  signal-present  hypothesis, 

logf 


”  =  D 

=  [9  (a^)]  • 


(87) 


For  the  Gaussian  shift  in  mean  signal  where  the  LLR  is  used  as  a  nonlinearity,  Broder’s 
asymptotic  performance  measure  is 


77  = 


2(72 


(88) 


Similarly,  such  asymptotic  relationships  can  be  derived  for  the  Page  test  with  nuisance 
parameter  estimation  for  the  Gaussian  shift  in  mean  signal.  First,  observe  that  Sk  of  equa¬ 
tion  (75)  may  be  approximated  by  the  first  term  in  the  infinite  summation  when  h  is  large 
enough, 

1  1  1 


5k  = 


(1  +  44)^’"*  {1+6 


+  ■  • 
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•# 
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r+l 


(89) 


Use  of  this  approximation  in  equations  (77)  and  (84)  results  in  the  large  threshold  approxi¬ 
mations 


Dw 


h-1 
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h 
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(90) 


and 


(91) 


where  rj  is  as  in  equation  (88)  and  it  has  been  assumed  that  J  is  also  somewhat  large.  The 
linear  relationship  between  h  and  Dw  is  seen  from  equation  (90).  However,  the  exponential 
relationship  that  exists  in  the  standard  Page  test  for  the  Gaussian  shift  in  mean  signal 
becomes  the  power  law  relationship  described  by  equation  (91)  when  the  nuisance  parameters 
are  estimated.  As  expected,  the  power  law  relationship  turns  into  the  exponential  one  of 
equation  (85)  when  the  amount  of  auxiliary  data  increases  to  infinity  (J  — >  oo), 

fw  (92) 

V 


This  may  be  seen  by  noting  that  [g  (x)]  =  r]  and  te^  =  1  for  the  Gaussian  shift  in 

mean  signal  case  using  the  LLR  and  by  applying  the  identity 

Urn  fl-^-y"  =  e“^  (93) 

\  xj 
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CONCLUSIONS 

This  report  has  presented  a  method  for  dealing  with  unknown  noise  or  nuisance  parameters 
in  the  Page  test  for  the  detection  of  the  onset  of  a  signal.  The  algorithm  is  a  generalization 
of  the  standard  Page  test  where  the  nuisance  parameters  are  estimated  using  data  previous 
to  a  most  recent  reset — data  declared  to  be  signal-free  by  the  Page  test.  This  also  makes 
analysis  feasible  due  to  the  independence  of  these  data  and  the  data  after  the  most  recent 
reset.  The  ASNs  of  the  Page  test  with  nuisance  parameter  estimation  were  derived  using 
both  Wald-  and  Siegmund-based  approximations.  It  was  shown  that  as  the  estimation  of  the 
nuisance  parameter  becomes  perfect  (i.e.,  infinite  auxiliary  data),  the  ASN  approximations 
of  the  Page  test  with  nuisance  parameter  estimation  simplify  to  the  ASN  approximations  of 
the  standard  Page  test. 

A  Gaussian  shift  in  mean  signal  with  an  unknown  variance  was  considered  as  an  example. 
Closed  forms  for  the  Wald-  and  Siegmund-based  approximations  to  the  average  time  before 
detection  and  between  false  alarms  for  the  Page  test  with  nuisance  parameter  estimation  were 
derived.  The  validity  of  the  approximations  was  verified  by  comparison  to  simulation  results 
where  it  was  observed  that,  as  in  the  standard  Page  test,  the  Siegmund-based  approximations 
provide  higher  accuracy.  The  Siegmund-based  approximations  were  used  to  investigate  the 
loss  in  performance  incurred  by  having  to  estimate  the  unknown  variance;  specifically,  the 
false  alarm  performance  (f)  as  a  function  of  threshold  and  the  detection  performance  (D) 
as  a  function  of  SNR. 

The  asymptotic,  in  the  sense  of  a  large  threshold,  relationship  between  the  threshold  and 
the  average  time  before  detection  for  the  Page  test  with  nuisance  parameter  estimation  was 
found  to  be  linear,  as  with  the  standard  Page  test.  However,  the  asymptotic  relationship 
between  the  threshold  and  the  average  time  between  false  alarms  was  found  to  follow  a 
power  law  that  approximates  the  exponential  relationship  of  the  standard  Page  test  when 
the  amount  of  auxiliary  data  used  to  estimate  the  nuisance  parameters  becomes  large. 
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